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The values of charged pion polarizabilities obtained in the framework of chiral per¬ 
turbation theory at the level of two-loop accuracy are compared with the experimental 
result recently reported by COMPASS Collaboration. It is found that the calculated 
value for the dipole polarizabilities (« — /3).„.± = (5.7± 1.0) X 10“^ fm® fits quite well the 
experimental result (a — /3),r± = l-2stat ± 1.4syst) x 10“^ fm®. 
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1. Introduction 

Almost ten years ago we have evaluated t he am plitude for 77 —>■ tttt in the framework 
of chiral perturbation theory (ChPT ) DEEl at two-loop order. It was done for 
neutral pions SI as well as for charged piones The obtained results were foun d to 
be in agreement with the only previous calculation performed at this accuracy 
provided that the sam e set of low-energy constants (LECs) is used. With updated 
LECs at order found ^ for the dipole polarizabilities the values 

[a - P)^± = (5.7 ± 1.0) X 10"^ W, 

{a -b /3)^± = 0.16 X 10”^ fm^ (1) 

At that time the MAMI Collaboration has reported the experimental result 

{a - /3)^± = (11.6 ± 1.5,tat ± ^-Osyst ± 0.5„,od) X 10“^ fm^ . (2) 

The index “mod” denotes the uncertainty generated by the theoretical models used 
to analyze the data. The ChPT calculation was clearly in conflict with the MAMI 
result, see also ^ for a recent discussion. 

The COMPASS collaboration at CERN has recently investigated pion Compton 
scattering by using the Primakoff effect. The pion polarizability has been determined 
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to be^ 


(a - = (4.0 ± 1.2,tat ± 1.4syst) X lO’^ fm^ (3) 

under the assumption (a+/3),n.± = 0. This result is in agreement with the expectation 
from chiral perturbation theory. 

The concept of the polarizability ofmolecules, atoms and nuclei was applied 
for the first time to hadrons in Refs. I 13 | 14 | 15 | gy using the general properties 
of quantum field theory it was shown that an expansion of the Compton scatter¬ 
ing amplitude for hadrons with spin one half in small photon energy up to the 
second order contains two structure parameters called the electric and magnetic 
hadron pol arizabilities. The classical sum rule for these quantities has been derived 
in Ref. Further theoretical investigation of the pion polarizabilities has been 
pursued since the early 1970s. In the current algebra + PCAC approach of Ter¬ 
ent’ev the fundamental low-energy theorem has been proven which allows one 
to relate the pion polarizability to the ratio 7 = /iyi( 0 )//iy( 0 ) of the vector and 
axial form factors in radiative pion decay tt —>■ ez/ 7 . By using recent precise mea¬ 
surements of the pion weak form factors by the PIBETA collaborationLl^one finds 
a^± = = 2.78(10) X 10-4 fm^ 

There were many calculations of the pio n polarizabilities by employing various 
models: the linear tr-model with quarks the chiral quark model 1201 ^ the super¬ 
conductor quark model some chiral models and so on. 

Almost all of them except Terent’ev approach predicted a value of the electric 
polarizability within the range 

4.0 X 10-4 fm^ < a^± < 6.0 x 10-4 fm^ (4) 

which we call large-valued results. We note that models not based on a chiral La- 
grangian, i.e., dispersion relations and hnite-e nergy sum rules, also obtained the 
polarizability within this range of values l ^3 | 24 | p ion and kaon polarizabilities 
have been calculated in the quark confinement model ^23 which the emphasis is 
placed on quark conhnement and the composite nature of hadrons. It was found for 
charged pions ^ 3.6 x 10-4 fm^ which is smaller than the large-valued results 
but slightly larger than Terent’ev’s prediction. 

The first correct calculation of the cross section 77 —>■ tttt within chiral perturba¬ 
tion theo ry t o next-to-leading order (one-loop accuracy) was performed in^^SI yy^s 
shown inl23^;Pa^; chiral symmetry relates the low-energy constants (LECs) appear¬ 
ing in the 77 —> TTTr-amplitude with the axial form factor /lyi(O). Thus it was shown 
explicitly that Terent’ev’s low-energy theorem follows from one-loop calculation of 
the 77 —>■ TTTT process within chiral perturbation theory. 

Note that the axial form factor 1 ia( 0) can be expressed through the dispersion 

jooi 

integral of the difference of the vector and axial spectral densities By using this 
sum rule the pion polarizability was estimated in ^23 and found to be in perfect 
agreement with chiral perturbation theory. 










September 18, 2015 0:43 WSPC/INSTRUCTION FILE HS-15.Ivanov 


Pion polarizabilities: Theory vs Experiment 3 


An actual two-loop ChPT calculation of the 77 —>■ tttt amplitude was done in 
(neutral pions) and® (charged pions). Because the effective Lagrangian at order 
p® was not available at that time, the ultraviolet divergences were evaluated in the 
MS scheme, then dropped and replaced with a corresponding polynomial in the 
external momenta. The three new counterterms which enter at this order in the 
low-energy expansion were estimated with resonance saturation. Whereas such a 
procedure is legitimate from a technical point of view, it does not make use of the 
full information provided by chiral symmetry. 

Later on, considerable progress has been made in this field, both in theo ry and 
experiment. As for theory, the Lagrangian at order has been constructed 

1001 

and its divergence structure has been determined This provides an important 
check on the above calculations: adding the counterterm contributions from the 
p® Lagrangian to the MS amplitude evaluated in ® and in ® must provide a scale 
independent result. Also in the theory, improved techniques to eval uate the two- 
loop diagrams that occur in these amplitudes have been developed!^. The updated 
calculation of the 77 —>■ tttt amplitude to two loops was then performed in® (neutral 
pions) and ® (charged pions). The final results for the pion polarizabities were 
presented in a rather compact algebraic form. By using updated values for the 
LECs one obtains the values of the pion polarizabilities given in Eq. O- 

rehensive review of the modern status of this field maybe found in 



Finally, one has to mention that research on pion polarizabities usin 



ulation is currently conducted by several groups, see, for instance. Refs. 

2. Definition of pion polarizabilities 

The electric (an) and magnetic {13h) polarizabilities characterize the response of 
hadron to two-photon interactions. These quantities must be considered as funda¬ 
mental as the electromagnetic mean square radii, static magnetic moments, etc. 
They are defined by the expansion of the Compton scattering amplitude in small 
photon momenta and energies. Since our interest here is the pion polarizabilities, 
we plot in Fig. [T] the diagram describing the Compton scattering by charged pion. 
Expanding the Compton scattering amplitude in small photon momenta and ener- 


Qi? Cl 


Q 2 , 62 


Pi 



P2 


Transverse gauge ej = (0, Ei) 
LAB frame pi = (M^^, 0) 


Fig. 1. A schematic illustration of Compton scattering by a pion. 
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-2e^ei • 6*2 

^ ^ > 

Born term 

+ 8 7rM.,r |a,7 wiW 2 ei ■ e 2 + (ei x qi) • (?2 x g 2 )| + • ■ • (5) 

el-mag polarizabilities 

It is convenient to use the linear combinations of the electric and magnetic polariz¬ 
abilities: (a — P)Tr and (a-I-/?),„• which are obtained from the helicity flip and helicity 
non-flip amplitudes, respectively. 

As follows from the definition, the dipole pion polarizabilities are proportional 
to a,n.(/3,n.) ^ ~ 4 X 10“'*fm^ where the hadronic scale A ^ dTrE^. ^ 1 GeV 

was used. Then a natural choice of units for the polarizabilities is 10“*fm^. 

3. Effective Lagrangian 

We consider an effective Lagrangian of QCD with two flavors in the isospin sym¬ 
metry limit m„ = md = m. At next-to-next-to-leading order (NNLO), one has 

El 

•^eff = ^2 + + -^6 ■ (6) 

The subscripts refer to the chiral order. The expression for £2 is 
/?2 

£2 = — + M^iJJ + U^)) , 

D^U = d^U -i{QU -UQ)A^, Q = |diag(l,-l), (7) 

where e is the electric charge, and denotes the electromagnetic held. The quantity 
F denotes the pion decay constant in the chiral limit, and is the leading term in 
the quark mass expansion of the pion (mass)^, = M^(l -|- 0(fh)). Further, the 
brackets (...) denote a trace in havor space. In Eq. 0, we have retained only the 
terms relevant for the present application, i.e., we have dropped additional external 
helds. We choose the unitary 2x2 matrix U in the form 

U = a + in/F, a2 + ^ = l2x2, (8) 

The Lagrangian at NLO has the structure^ 

7 3 . 

Ci = Y,hK, + Y,h^K, = ^{D^UD^^U^f + ■■■, ( 9 ) 

i=\ i=l 

where li ,hi deno te low-energy couplings, not hxed by chiral symmetry. At NNLO, 
one hastaM^ 

57 

£5 = CiPi . 
i=l 


gies, one hnds 

£y7r+—>7-7r'+ 


( 10 ) 
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As was shown in Ref. SSI the number of operators Pi can be reduced by at least one 
from 57 to 56. For the e xplicit exp ressions of the polynomials Ki,Ki and Pi, we 
refer the reader to Refs. vertices relevant for 77 —t tt+tt" involve 

li,... ,le from £4 and several c^’s from Cq, see below. 

The couplings h and Ci absorb the divergences at order and p®, respectively, 

^ -7f’( 7 !^ +7|^A.<^)) a} , 

The physical couplings are Z[(p,4) and c^(p,4), denoted by l ^, c[ in the following. 
The coefficients 7 ^ are given in^l^ and 7 !^’^’^^ are tabulated in SSI, yvb shall use the 
scale independent quantities li introduced in^l^ 

= + ( 12 ) 

where the chiral logarithm is I = ln(M^/p^). We shall use SI 

Ii = -0.4 ± 0.6 , I 2 = 4.3 ± 0.1, ^3 = 2.9 ± 2.4 , U = 4.4 ± 0.2 , (13) 

and 

[a = k - k = 3.0 ± 0.3 (14) 

obtained from radiative pion decay to two loop accuracy SSH. 

The constants occur in the combinations 

a{ = — 409671 ^ (6 Cg + C 29 — Cgg — 3 C 34 + C 35 + 2 c^q — 4 C 47 + Cgg) , 

al = 2567r4(8c^g -Sc^o +0^2 “ 2 + 4 c(j4 + 8 c^g - 4 c^^) , 

6 ’' = -1287r4(cSi+cS2-2c53-4c^4) . 

As follows from the resonance exchange model SI 

(^al,a 2 ,b^^ = ^ — 3.2,0.7, 0.4^ [ present work ] . (15) 

The values of these constants were obtained in the ENJL model SH {a\,a 2 ,b'^) = 
(—8.7, 5.9, 0.38) One can see that only P agrees in the two approaches. We shall 
use P = 0.4 ± 0.4. The combinations (a ± i3)t^± are independent of 02 and are 
determined precisely by the chiral expansion to two loops, once is fixed. We will 
then simply display this quantity as a function of - the result turns out to be 
rather independent of its exact value. 
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4. Evaluation of the diagrams 

The lowest-order contributions to the scattering amplitude are described by tree- 
and one-loop diagrams. These contributions were calculated in l2Sl xhe two-loop 
diagrams are displayed in Figs. [21 [4] and jS] The two-loop diagrams in Fig. [2] may 
be generated according to the scheme indicated in Fig. |3l where the filled in blob 
denotes the d-dimensional elastic TTTr-scattering amplitude at one-loop accuracy, 
with two pions off-shell. 

The diagrams shown in Fig. [4] may be reduced to tree-diagrams by using Ward 
identities. They sum up to the expression 





qiT{‘^P2 



Z^R{t) 


+ crossed 


(16) 


where Z^^ is the pion renormalization constant. The function R{t) starts at order 
l/F^ and can be obtained from the full pion propagator. 

Two further diagrams are displayed in Fig. |SJ The first one - called “acnode” 
in the literature - may again be evaluated by use of a dispersion relation, see SI . 
The second one is trivial to evaluate, because it is a product of one-loop diagrams. 



Fig. 2. A set of two-loop diagrams generated by C2 and one-loop diagrams generated 
by £4. 
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The remaining diagrams at order p® are shown in Fig. |6l 

The evaluatio n o f the diagrams was done in the manner described in EI331 i^y 
invoking FORM S31. in particular, we have verihed that the counterterms from 
the Lagrangian Cq remove all ultraviolet divergences, which is a very non-trivial 
check on our calculation. Furthermore, we have checked that the (ultra-violet finite) 
amplitude so obtained is scale independent. 



Fig. 3. Construction scheme for the diagrams in Fig. [2l 


9i Pi 9i Pi 



Fig. 4. A class of one-particle reducible diagrams. The filled in circles summarize self¬ 
energy and vertex corrections. 



Fig. 5. Acnode and butterfly diagrams. 
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( 1 ) ( 2 ) ( 3 ) ( 4 ) 


Fig. 6. The remaining diagrams at order p®: one-loop graphs generated by £4, and coun¬ 
terterm contributions from Tg. 


5. Chiral expansion for pion polarizabilities 

Using the same notation as in^ we find for the dipole polarizabilities 


{a ± / 3 )^+ = 


a 


16 7r2 F2 


Cl± 


16 7r2 F2 


0(M4) , 


where 


2 - 

ci+ = 0, Cl- = -Ia , 


4 r, 

/, 1- 3-\ 

53 , 

1 - 

3 r 

91 . 1 


( ^ 7: 7: ^2 ) 

- + 

n 1 4” 

nh + 

+ A_(. ? 

9 1 

V 2 2 J 

24 

2 

2 

72 


65 


di_=a\ + 8b^--U[h-l2 + lA-^]-::li-^l2 + ^h- I Ah 


12 


1 r 1 r 


187 ^ 

--h A- 

108 


with 


A+ = 


8105 135 


A 41 53 2 

A- =-TT^ . 

432 64 


(17) 


(18) 


(19) 


576 64 

It would be interesting to numerically compare the values of A± given by Eq. (HU) 
with those obtained in Refs. One has 


^ _ f —6.75 our ^ _ f —8.08 our 

+ “ \ -8.69 Burgi " “ (-8.73 Burgi 

The results for the polarizabilities evaluated with the central values for the 
LECs in Eqs. (fT3l) - (IT5]) are shown in Table [T] The uncertainty in the prediction for 


Table 1. Central values of polarizabilities in 
units of 



to one loop 

to two-loops 

+ 

1 

6.0 

5.7 

{a + d)n+ 

0 

0.16 


the polarizability has two sources. First, the low-energy constants are not known 
precisely. Second, we are dealing here with an expansion in powers of the momenta 
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and of the quark masses up to and including terms of order p®. The discussion of 
estimating uncertainties may be found in our paper It was shown that the value 
for the dipole polarizability (a — /3),r± is rather reliable - there is no sign of any 
large, uncontrolled correction to the two-loop result. The maximum deviation 1.0 
from the central value 5.7 has been used as the final theoretical uncertainty for the 
dipole polarizability: 

(a-/3)^± = (5.7±1.0) X lO-'^fm®. (20) 

The chiral expansion for the combination (a -I- I3)t^± starts out at order p® so we 
have determined only its leading order term. 


6. Experimental information 

There are three types of experiments aiming to measure the pion polarizabilities: 

• The scattering of high energy pions off the Coulomb field of heavy nucleus 
using the Primakoff effect 

• Radiative pion photoproduction from the proton 

• Pion pair production in photon-photon collisions 

Schematically, they are shown in Fig. [T] The possibility to measure the pion 
polarizability via the Primakoff reaction was proposed in the early 1980s in EH. 
The measurement of the pion-photon Compton scattering amplitude by using the 
Primakoff effect was performed in an experiment at Serpukhov EH ^ but the small 
data sample led to only an imprecise value for the polarizability of = (6.8 ± 
1.4stat±l-2syst) X 10“"^ fm®. Low statistics made it difficult to evaluate the systematic 
uncertainty. 

COMPASS has now achieved a modern Primakoff experiment, using a 190 GeV 
pion beam from the Super Proton Synchrotron at CERN directed at a nickel target. 
It is important that COMPASS was also able to use a muon, which is point-like par¬ 
ticle, to calibrate the experiment. The Compton 7r“7 —>■ 7r“7 scattering is extracted 
from the reaction 7r“ Ni —^ Ni by selecting events from the Coulomb peak at 

small momentum transfer < 0.0015 GeV^. From the analysis of a sample of 



(a) (b) (c) 


Fig. 7. Schematic images of the available experiments. 
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63,000 events, the collaboration obtained a value of the pion electric polarizability 
of Ott = ( 2.0 ± O.Gstat ± O.Tsyst) X 10 “^ fm^ under assumption {a + /3)^ = 0 . 

The cross section for the radiat ive pion photoproduction yp —> yTr+n has been 
measured at the Lebedev Institute By using an extrapolation to the pion pole 
in the unphysical region the value of the electric polarizability was obtained = 
—/3^+ = (20 ± 12) X 10“^ fm^. Similar experiment was performed at the Mainz 
Microtron MAMI^J^but the pion polarizability has been extracted by a comparison 
of the data with the predictions of two different models yielding the value {a — 
/?),+ = (11.6 ± 1.5stat ± S.Osyst ± 0.5„,od) X 10-4 fm^. 

Another possibility to obtain the value for the pion polarizability is to extrap¬ 
olate the data from the pion pair production in photon-photon collisions 77 —> tttt 
to the region of the Compton scattering threshold by using crossing symmetry and 
analyticity. Normally, the procedure involves the construction of the dispersion re¬ 
lations with one or two subtractions. The most recent analysis preformed in ^231 
produced the value (a — /3)^± = Id.O^i^g x IQ-^fm^ which is close to the MAMI 
data. There are plenty of previous studies in this direction which give quite a broad 
region for the value of the pion polarizability. The available experimental informa¬ 
tion is shown in Table [2] 


Table 2. Experimental informatio n on (a — I3) ^±, in units of 10 "^fm®. We indi¬ 
cate the reaction and data used. In EZl ^ HU and El was determined, using as a 
constraint o.,r = —P-n- To obtain (o — , we multiplied the results by a factor of 

2 . 


Experiments 

(« - I3)^± 

'K~ Z ^ Z Serpukhov (1983)1^ 

15.6 ± 6.4atat ± 4.4syst 

IP 77r"*'n Lebedev Inst. (1984) 

40 ±24 

D. Babusci et al. (1992) 

77 TT+TT- PLUTO ESI 

DM 1 HSl 

DM 2 ESl 

mark II ED 

38.2 ±9.6 ±11.4 

34.4 ±9.2 

52.6 ±14.8 

4.4 ±3.2 

J.F. Donoghue, B. Holstein (1993) 

77 TT+TT- MARK IlEH 

5.4 

A. Kaloshin, V. Serebryakov (1994) 

77 -)• TT+TT- MARK II ES 

Crystal Ball Coll. EH 

5.25 ±0.95 

IP 77r"*'n Mainz (2005) 

11.6 ± 1.5stat ± S.Osyst ± 0.5 mod 

L. Fil’kov, V. Kashevarov (2005) 1^ 

77 ^ TT+TT- MARK II ED, 

TPC/ 27 E 5 I, CELLO ESl^ 

VENUS E3, ALEPH ES, BELLE EH 

IS.Otjg 

TT-Z ->■ 7TT-Z COMPASS (2015)Ej 

4.0 ± 1.2stat ± 1.4syst 
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